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Top-heavy IMFs in ultra-compact dwarf galaxies? 
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Abstract. Ultra compact dwarf galaxies (UCDs) are dense stellar systems at the 
border between massive star-clusters and small galaxies. Their on average high optical 
mass-to-light (M/L) ratio cannot be explained by stellar populations with the canonical 
stellar initial mass function (IMF), while it is doubtful that non-baryonic dark matter can 
accumulate enough on the scales of UCDs for influencing their dynamics significantly. 
UCDs in the Virgo galaxy cluster apparently also have an over-abundance of neutron 
stars, strongly suggesting a top-heavy IMF, which would explain both findings. This is 
because a top-heavy IMF can provide the unseen mass through an abundance of stellar 
remnants. The suggested variation of the IMF can be understood if UCDs represent 
a case of rapid star-formation in an extremely dense environment. While top-heavy 
IMFs imply a much heavier mass-loss shortly after the formation of a stellar system, 
this process does not necessarily dissolve the UCDs. Their formation with a top-heavy 
IMF would therefore not contradict their existence. 



1. Introduction 

Ultra-compact dwarf galaxies (UCDs) are stellar systems with V-band luminosities be- 
tween 10 6 and some 10 7 Solar luminosities, but half-light radii of only about 50 pc 
or less. Most confirmed UCDs are found in nearby galaxies with rich globular clus- 
ter systems, such as NGC 5128 (Centaurus A), or nearby galaxy clusters, such as the 
Fornax cluster and the Virgo cluster. Even with the most powerful telescopes available 
at present, UCDs are barely resolved at those distances. Their true nature has there- 
fore only been revealed quite recently by estimating their distance from their redshift 
(Hil ker et al.lll999l : lDrinkwater et al.ll2000l) . Before that, UCDs were either interpreted 



as background galaxies or foreground stars, because of their appearance as bright point 
sources in most observations. Colours and in some cases also line indices suggest 
intermediate- to high ages for UCDs. 

Due to the distance and the compactness of UCDs, it is difficult to get information 
on their internal stucture. The observed density profiles of UCDs carry a strong imprint 
of the point-spread function and the spectrum of a UCD is obtained in observations that 
cover a good fraction of the whole UCD. Estimates of their central velocity dispersions, 
their global velocity dispersions, their core radii or their half- light radii therefore require 
involved modelling of their density profiles, as done by, e.g. iHilker et"aT1 d2007l) . 



A comprehensive list of globular clusters (GCs) and UCDs for which such mod- 
ell ing has been done, l eading also to estimates for their dynamical masses, M, is found 



m 



Mieske et all (120081) . The half-light radii (r^) and dynamical mass-to-light (M/L) 
ratios of the objects in that list are shown in Fig. [J According to this figure, UCDs are 
indistinguishable from GCs at a a luminosity of « 10 6 L©, corresponding to a dynamical 
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Figure 1 . Half-light radii (left panel) and dynamical M/L ratios (right pan el) of 
GCs and UCDs, using the compilation of such objects by Mies ke et alj y008). The 
dashed vertical lines at a dynamical mass, M, of 2 x 10 6 M o mark a definition that 
seperates GCs from UCDs, because objects with masses higher than 2x 10 6 M show 
a clear tendency to higher and M/L ratios, even though GCs are indistinguishable 
from UCDs at M m 2 x 10 6 M Q , Note that also the two-body relaxation times become 
longer than a Hubble time M su 2 x 10 6 M , i.e. GCs are relaxed stellar systems and 
UCDs are not. 



mass of 2 x 10 6 M Q (see also Mi eske et al.ll2008). This has led some authors to consider 
UCDs as the most massive GCs (e.g. Mieske et al.ll2002r) or to re-interpret oj Cen, the 
most massive GC of the Milky Way, as a low-mass UCD. However, the average rh) 
and average dynamical (M/L) ratios of objects with M > 2 x 10 6 M Q increase system- 
atically with M, whereas this is not the case for objects with M < 2 x 10 6 M Q . This 
motivates to define the objects with M > 2 x 10 6 M Q as UCDs, in order to distinguish 
them from 'classical' GCs. 

The origin of UCDs is still a debated issue. However, numerical experiments pre- 
dict that the young and massive star-cluster complexes observed in interacting galaxies 
evolv e into objects very much like a UCD. Th e time required for this is only a few 100 
Myr (Kroupa 1998: iFellhauer & Kroupall2002l) . Supporting evidence for this notion is 
provided by the discovery of W3, an object which was probably created by the forma- 
tion of the merger-remn ant galaxy NGC 7252 and which is a UCD in ter ms of its mass 
and its half-light radius dMaraston et al.ll2004l : IFellhauer & Kroupall2005h . 



2. The M/L ratios of UCDs 

Even if GCs and UCDs have the same origin, their differences in M and r\ imply that 
they have evolved differently, which has implications on the difference of the average 
M/L ratio between GCs and UCDs. The key to understanding this is that each stellar 
system is subject to two kinds of evolution: Stellar evolution and dynamical evolution. 
Stellar evolution implies that the brightest stars evolve first and thereby become dark 
remnants, leading to an increase of the M/L ratio of a stellar system, unless they are 
ejected. Dynamical evolution, on the other hand, tends to lower the M/L ratio of a 



Top-heavy IMFs in UCDs 



3 



stellar syst em, as it leads to a preferent ial loss of low-mass stars, which have a high 
M/L ratio dBaumgardt & Makinol |2003). It is evident that the time-scale on which 
stellar evolution changes the properties of a stellar system is independent of M and 
rh of the stellar system, but this is not the case for dynamical evolution. A measure 
for whether a stellar system is affected by dynamical evolution is its m edian two-bod y 
relaxation time, t r h, which can be estimated using equation (2-63) from lSpitzed (Il987h . 
If t r f, is longer than the age of the system, it can be considered unaffected by dynamical 
evolution. 

Using their present-day parameters, it turns out that GCs usually have ? r h very 
much b elow a Hubble time, where as UCDs have of the order of a Hubble time or 
longer (Dabringhausen et al.ll2008h . The M/L ratios of GCs are thereby determined 
by an interplay of their s tellar evolution as well as their dynamical evolution (see 
Kruijssen & Mieskell2009l) . UCDs are, in comparison, rather simple systems, as the 



evolution of their M/L ratio is essentially determined by stellar evolution. Therefore, 
their M/L ratio has only increased after their star formation came to an end. Thus, as 
the stellar populations of GCs and UCDs generally are of the same age (with the ex- 
ception of W3), the average M/L ratio of UCDs being higher than that of GCs appears 
natural. 

As UCDs are thus nearly unaffected by dynamical evolution, their present-day 
stellar population should be determined by their IMF and their star formation history. 
A good initial assumption for the IMF of UCDs is the canonical IMF, since this IMF 
is fou nd to be consistent with all currently resolved stellar populations (lKroupall200ll 



2008). However, the observed dynamical M/L ratios of a clear majority of UCDs ex- 



ceed the predictions for a s ingle-age, single-metallicit y s tellar population (S SP) with 
the canonical IMF (see e.g lDabringhausen et al.ll2008l or LMieske et all 120081) . This is 
even the case if the stellar populations of all UCDs are assumed to be 13 Gyr old SSPs, 
i.e. if all stars in UCDs were nearly as old as the Universe. Note that for a given IMF 
and metallicity, no stellar population can have a higher M/L ratio than this oldest possi- 
ble stellar population. Also note that UCDs may consist of multiple stellar populations, 
as is the case with u Cen, and that the sub-populations in them may have different 
metallicities, as well as different ages. However, the metallicity of a UCD is estimated 
from its integrated light, leading to some value that should be representative for all stars 
in the UCD. 



3. Top-heavy IMFs in UCDs? 

3.1. Evidence for a top-heavy IMF in UCDs from their M/L ratios 

The surprisingly high M/L-ratio of UCDs implies some unseen mass in them. This can- 
not be the non-baryonic dark matter that is often invoked to explain, e.g., the rotation 
curves of spiral galaxies or the M/L-ratio of galaxy clusters. This is bec ause structures 



of co l d dark matter h ave been shown to obey a universal density profile (INavarro et al. 



19971 : iBullock et al.ll2001r) . This density profile implies that only a very small amount 



of cold dark matter gathers within the quite compact UCDs and their M /L-ratios soul d 
therefore be indistinguishable from the ones of pure stellar populations (Murray 2009). 
The internal accelerations of UCDs are also above the transition between pure Newto- 
nian dynamics and MOND, which is a modification of the classical theory of gravitation 
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Figure 2. The high-mass IMF slopes of UCDs (cf. Eq.[T|l, as suggested by their 
M/L ratio as a function of their observed over their predicted M/L ratio, assuming 
the canonical IMF (left panel), and as a function of their dynamical mass (right 
panel). For the curve in the left panel, it is assumed that 20 per cent of the NSs are 
retained, which seems to be a reasonable assumption (see Section lXTT i. The average 
value for the M/L ratio of UCDs over the prediction for the canonical IMF (dotted 
vertical line) with its l-cr variance (shaded area) suggests that the IMF in UCDs is 
considerably top-heavy, even though the assumtion that they all are 13 Gyr old is 
very conservative. 



in the limit of weak gravitational fields (|Milgromlll983|) . Thus, also MOND cannot ex- 
plain the M/L-ratios of UCDs. 

Excluding those possiblities leads to the assumption of an IMF that is not univer- 
sal, but varying, as an explanation for the high M/L ratios of UCDs. T his notion of a 
varying IMF does not neccesarily contradict the results by iKroupal (l200lh . as these re- 
sults were obtained considering resolved stellar populations in the Milky Way; i.e. from 
stellar populations in our immediate neighbourhood. UCDs, however, may have formed 
in much stronger starbursts than the ones that have occured in the Local Group. These 
extreme initial conditions would then be what distinguishes UCDs from the open clus- 
ters and GCs in the Milky Way. It has indeed been shown t hat extremely hig h densities 
would lead to stellar collisions, which would alter the IMF (|Bonnell & Batell2002l) . 

A varying IMF can be formulated as 

£0ra*) = kkimZ"', (1) 

with 

a x = 1.3, 0.1 < — <0.5, 

M Q 

/77 . 

a 2 = 2.3, 0.5 < — ^ < 1, 

M Q 

a 3 e E, 1 < — < 150, 



M. 



© 



where «i» is the initial stellar mass, the factors k[ ensure that the IMF is continuous 
where the power changes and k is a normalisation constant (iKroupa| [2008). £(m t ) equals 
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Figure 3. Comparison of the encounter rate to the fraction of GCs and UCDs in 
the Virgo galaxy cluster containing an LMXB (left panel) and the high-mass IMF 
slope that can explain the abundance of LMXBs in UCDs (right panel). The x-axes 
in both panels have been re-scaled, so that they show the average dynamical mass 
corresponding to a certain luminosity instead of the luminosity itself. This allows 
easy comparisons between this figure and Figs. Q] and [2] The solid line in the left 
panel shows the encounter rate under the assumption that the IMF in GCs and UCDs 
is canonical. By assuming that the IMF in GCs is canonical, but given by Eq. Q] 
in UCDs, the encounter rate in UCDs can by changed into the one indicated by the 
dashed line. This requires a certain dependency of the high-mass IMF slope on the 
mass of the UCDs. This function is shown in the right panel. It is remarkably similar 
to the one shown in the right panel of Fig. [2] even though it was derived entirely 
differently. 



if m, < 0.1 Mo or m» > 150 Mp, where 150 M^ arguably is the upper mass-limit 
for stars (e.g. IWeidner & Kroupal 120041 : lOev & Clarke! 120051) . Eq.Q] is the canonical 
IMF for «3 = 2.3. For < 2.3, the IMF is top-heavy, implying more intermediate- 
mass stars and in particular more high-mass stars. In old stellar populations like the 
UCDs, these stars have evolved into white dwarfs and neutron stars (NSs) and thereby 
contribute to the mass of the total stellar population but hardly to its luminosity. We 
note that the IMF-variation proposed here is not the o nly one that can explain high 
M/L rat i os. Al ternative scenarios have been discussed by Mies ke & Kroupal d2008l) and 
iMurrayl (2009). However, the advantage of the top-heavy IMF used here is that it can 
explain several properties of UCDs at the same time. 

Given Eq. [TJ it is possible to translate some observed M/L ratio divided by the 
M/L-ratio predicted by a SSP-model into some value of 0-3 that corresponds to this 
ratio by providing the appropriate amount of dark matter in form of WDs and NSs (see 
left panel of Fig. |2]). Note that the fraction of NSs that actually remain bound to the 
UCD is an im portant issue when calc ulating a^. Given the large peculiar velocities of 
many pulsars dLyne & Lorimer|[l994 ). which are thought to be NSs, a lot of NSs may 
escape from the UCDs. On the other hand, it is observationally confirmed that GCs 
and especially UCDs do have a population NSs (visible as low-mass X-ray binaries, cf. 
Section [X2]>. 
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The right panel of Fig. Q] shows that there is a relation between the M/L of UCDs 
and their dynamical mass. Combining this information with the as implied by a certain 
M/L ratio (left panel of Fig. [2]) leads to a relation betwe en the average of UCDs as 
a function of their mass (see IDa bringha usen et al.ll2009l for details). This relation, as 
shown in fig. 5 in IDabringhausen et al.1 d2009h . is plotted in the right panel of Fig. |2] 

3.2. Evidence for a top-heavy IMF in UCDs from their low-mass X-ray binaries 

If UCDs have indeed formed with the IMF formulated in Eq. \T\ and the IMF is top- 
heavy (i.e. Qf3 < 2.3), they should have a large population of NSs. Thus, compared to 
the case with the canonical IMF, UCDs would also have a higher number of close binary 
systems of a low-mass star and a NS. In such systems, the NS can accrete matter from 
its companion star and thereby become a bright X-ray source; a so-called low-mass 
X-ray binary (LMXB). The NSs in LMXBs thereby become observable, even though 
they remain invisible in the optical. 

Observations show that the densest stellar systems also have the highest likelihood 
to harbour a bright LMXB. This suggests that they are formed by three-body encoun- 
ters; especially encounters between a star and a binary. This is because these en counters 
tighte n binary systems if the binary entering the interaction was already tight (Heggie 
Il975h . It thus appears that stellar dynamics is the most important and perhaps even the 
only important facto r for the formation of LMXBs (IVerbuntll2003l; ISivakoff et alj|2007l ; 



Peacock et aDl2010h . 



The number of encounters relevant for the creation of LMXBs, i.e. encounters that 
involve stars and NSs, can be written as 



r oc C , (2) 

0~ 



where « ns is the number density of NSs, n s i s the number d ensity of other stars, r c is the 
core radius and cr is the velocity dispersion (IVerbuntH2003r) . Rather than measuring the 
volume of the stellar system by r^, we use rjj. This is a quantity that is available for all 
GCs and UCDs plotted in Fig. Q] and arguably only changes the proportionality factor 
in Eq. |2] 

Interesting for the purpose here is to find the average encounter rate of GCs and 
UCDs as a function of their luminosity, T(L). This quantity can be compared to the 
fraction of GCs and UCDs in certain luminosity-bins that actually have an LMXB, P, 
which is an estimator for the probability to find an LMXB in a system in the according 
luminosity range (and a corresponding mass range). Observational data on which frac- 
ti on of GCs and UCDs in the Virgo-galaxy cluster have a bright LMXB can be found 
in Sivako ff et all ((2007). If the stellar mass function in GCs is the same as in UCDs, the 
shape of Y only depends on how the average of GCs and UCDs changes with their 
luminosity, which is a function that can be derived using the data plotted in Fig. Q] A 
non-changing stellar mass function for all GCs does indeed suggest that the observed 
rise of P matches the increase of T expected for them. However, P still rises for UCDs, 
whereas the calculated T drops steeply if they all had the same stellar mass function (cf. 
left panel of Fig. [3]). 

This leads to the alternative assumption that the mass function of UCDs does 
change with their luminosity. An IMF that becomes increasingly top-heavy with the 
luminosity of UCDs implies that n s in Eq. [2] changes with luminosity. Consequently, 
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also cr increases, because more NSs imply a higher total mass. Assuming the popula- 
tion of NSs and other stars in GCs is given according to the canonical IMF and that the 
IMF in UCDs is given by Eq. \T\ T in UCDs can be calculated such that it is not only 
proportional to P in GCs, but also in UCDs. This implies that that the high-mass slope 
of the IMF changes with their luminosity (mass), as shown in the right panel of Fig. |3] 
Comparing the right panel of Fig. [3] with the right panel of Fig. |2] shows a qualitative 
agreement of how the high-mass IMF slope of UCDs would have to change in order 
to explain their M/L ratios and how it would have to change in order to explain the 
abundance of LMXBs in them (Dabringhausen & Kroupa, in preparation). 



4. Implications for the inital conditions in UCDs 

The dynamical evolution of young stellar systems is strongly influenced by rapid mass- 
loss, especially if their IMF is top-heavy. This is because of the strong radiation of 
massive stars, which can quickly remove all the interstellar medium that is not used up 
in star formation. Massive stars evolving to SN can in principle replenish the interstellar 
medium, but provide at the same time so much energy that also this gas is more likely 
to leave the stellar system, even if it is as massive a UCD. Stellar evolution can thus 
lead to further mass-loss that becomes more rapid and more severe with the the top- 
heaviness of the IMF. This mass-loss would lead to a strong expansion or even to a 
dissolution of the stellar system. The extreme IMFs suggested for UCDs in Sections [3TT1 
and 13.21 therefore imply that UCDs must have b een extremely compact if they are to 



survive until today with the observed densities dDabringhausen et aDl2010l) . Assuming 
that their stellar populations formed on a time-scale of a few Myr (as implied by the 
formation scenario suggested in Section [T), the more massive UCDs may have had 
initial masses of some 10 8 M Q while having initial half-light radii of only a few parsec. 
A top-heavy IMF would then imply a population of the order of 10 6 O-stars with a 
total luminosity of the order of 10 11 L Q . If this is the case, UCDs should be bright 
point-like sources, identifiable by a spectrum that is consistent with a stellar population 
(Pflamm-Altenburg, Dabringhausen & Kroupa, in preparation). 

The high initial density implied by the initial mass and initial radius could explain 
why the IMF was top-heavy in the first place. At such densities, collisions between 
proto-stars become co mmon, which may cause a deviation from the canonical IMF 
(iBonnell & Batell2002h . 



5. Conclusion 

The observed dynamical M/L ratios of UCDs and the number of LMXBs in them can 
be used to constrain their IMF. Their high MIL ratios are a strong argument for a non- 
canonical IMF in UCDs, but the MIL ratios alone leave a number of possibilities on how 
the IMF could vary with the mass of the UCDs. However, UCDs apparently also have a 
richer population of NSs than GCs. This cannot be explained by simply assuming that 
more NSs remain bound to UCDs, while GCs and UCDs both formed with the canonical 
IMF, because in this case the M/L ratios of UCDs should not tend to exceed the value 
predicted for a canonical IMF. Thus the number of LMXBs in UCDs and the M/L 
ratios of UCDs in combination imply that the IMF of UCDs was top-heavy, especially 
for very high-mass stars. This implies that UCDs would have formed compact and 
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with very high densities. Note that the number of LMXBs also give a constraint on 
the low-mass IMF in UCDs, because low-mass stars have to be rather frequent as well. 
Otherwise, LMXBs in UCDs would be rare for a lack of low-mass companions for the 
NSs (cf. EqO. 
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